Forging the path to precision medicine in Qatar: a public health perspective on pharmacogenomics initiatives.

Pharmacogenomics (PGx) is an important component of precision medicine that promises tailored treatment approaches based on an individual's genetic information. Exploring the initiatives in research that help to integrate PGx test into clinical setting, identifying the potential barriers and challenges as well as planning the future directions, are all important for fruitful PGx implementation in any population. Qatar serves as an exemplar case study for the Middle East, having a small native population compared to a diverse immigrant population, advanced healthcare system, national genome program, and several educational initiatives on PGx and precision medicine. This paper attempts to outline the current state of PGx research and implementation in Qatar within the global context, emphasizing ongoing initiatives and educational efforts. The inclusion of PGx in university curricula and healthcare provider training, alongside precision medicine conferences, showcase Qatar's commitment to advancing this field. However, challenges persist, including the requirement for population specific implementation strategies, complex genetic data interpretation, lack of standardization, and limited awareness. The review suggests policy development for future directions in continued research investment, conducting clinical trials for the feasibility of PGx implementation, ethical considerations, technological advancements, and global collaborations to overcome these barriers.

Family History of Twinning and Fertility Traits in Nigerian Mothers of Dizygotic Twins.

Familial twinning and fertility traits were investigated in Nigerian mothers of dizygotic (DZ) twins (MoDZT; N = 972) and controls (N = 525) who responded to our person-to-person interview, which included questions on pregnancy history and family history of DZ twinning. Controls were defined as women who are not twins themselves and do not have twins in their first-degree relatives. Over 95% of the participants were Yoruba. We found that Nigerian MoDZT had an average of 4.0 (±2.6) pairs of twins among their relatives, and of these, the prevalence of DZ twins was significantly higher than that of monozygotic (MZ) twins (45.9% vs. 25.8%). Controls had an average of 0.5 (±0.4) pairs, and over 95% of the controls had no twins in their relatives. These results suggest genetic influences on DZ twinning in Nigerians. MoDZT were significantly younger in their mean age at first child, and had higher parity than controls, suggesting increased fertility in MoDZT. As compared to mothers with a single set of twins, mothers (N = 130) with multiple sets had significantly more twins among their relatives (5.4 pairs vs. 3.7 pairs) and had their first twins at a younger age (28.4 vs. 30.7 years), indicating that mothers with multiple sets of twins might have higher genetic propensity for twinning associated with earlier age at twin pregnancy. Our findings argue for genomewide association studies for DZ twinning in Nigerians, and may help to develop intervention strategies to overcome infertility/subfertility problems.

Genome-wide association study meta-analysis of dizygotic twinning illuminates genetic regulation of female fecundity.

STUDY QUESTION: Which genetic factors regulate female propensity for giving birth to spontaneous dizygotic (DZ) twins? SUMMARY ANSWER: We identified four new loci, GNRH1, FSHR, ZFPM1, and IPO8, in addition to previously identified loci, FSHB and SMAD3. WHAT IS KNOWN ALREADY: The propensity to give birth to DZ twins runs in families. Earlier, we reported that FSHB and SMAD3 as associated with DZ twinning and female fertility measures. STUDY DESIGN, SIZE, DURATION: We conducted a genome-wide association meta-analysis (GWAMA) of mothers of spontaneous dizygotic (DZ) twins (8265 cases, 264 567 controls) and of independent DZ twin offspring (26 252 cases, 417 433 controls). PARTICIPANTS/MATERIALS, SETTING, METHODS: Over 700 000 mothers of DZ twins, twin individuals and singletons from large cohorts in Australia/New Zealand, Europe, and the USA were carefully screened to exclude twins born after use of ARTs. Genetic association analyses by cohort were followed by meta-analysis, phenome wide association studies (PheWAS), in silico and in vivo annotations, and Zebrafish functional validation. MAIN RESULTS AND THE ROLE OF CHANCE: This study enlarges the sample size considerably from previous efforts, finding four genome-wide significant loci, including two novel signals and a further two novel genes that are implicated by gene level enrichment analyses. The novel loci, GNRH1 and FSHR, have well-established roles in female reproduction whereas ZFPM1 and IPO8 have not previously been implicated in female fertility. We found significant genetic correlations with multiple aspects of female reproduction and body size as well as evidence for significant selection against DZ twinning during human evolution. The 26 top single nucleotide polymorphisms (SNPs) from our GWAMA in European-origin participants weakly predicted the crude twinning rates in 47 non-European populations (r = 0.23 between risk score and population prevalence, s.e. 0.11, 1-tail P = 0.058) indicating that genome-wide association studies (GWAS) are needed in African and Asian populations to explore the causes of their respectively high and low DZ twinning rates. In vivo functional tests in zebrafish for IPO8 validated its essential role in female, but not male, fertility. In most regions, risk SNPs linked to known expression quantitative trait loci (eQTLs). Top SNPs were associated with in vivo reproductive hormone levels with the top pathways including hormone ligand binding receptors and the ovulation cycle. LARGE SCALE DATA: The full DZT GWAS summary statistics will made available after publication through the GWAS catalog (https://www.ebi.ac.uk/gwas/). LIMITATIONS, REASONS FOR CAUTION: Our study only included European ancestry cohorts. Inclusion of data from Africa (with the highest twining rate) and Asia (with the lowest rate) would illuminate further the biology of twinning and female fertility. WIDER IMPLICATIONS OF THE FINDINGS: About one in 40 babies born in the world is a twin and there is much speculation on why twinning runs in families. We hope our results will inform investigations of ovarian response in new and existing ARTs and the causes of female infertility. STUDY FUNDING/COMPETING INTEREST(S): Support for the Netherlands Twin Register came from the Netherlands Organization for Scientific Research (NWO) and The Netherlands Organization for Health Research and Development (ZonMW) grants, 904-61-193, 480-04-004, 400-05-717, Addiction-31160008, 911-09-032, Biobanking and Biomolecular Resources Research Infrastructure (BBMRI.NL, 184.021.007), Royal Netherlands Academy of Science Professor Award (PAH/6635) to DIB, European Research Council (ERC-230374), Rutgers University Cell and DNA Repository (NIMH U24 MH068457-06), the Avera Institute, Sioux Falls, South Dakota (USA) and the National Institutes of Health (NIH R01 HD042157-01A1) and the Genetic Association Information Network (GAIN) of the Foundation for the National Institutes of Health and Grand Opportunity grants 1RC2 MH089951. The QIMR Berghofer Medical Research Institute (QIMR) study was supported by grants from the National Health and Medical Research Council (NHMRC) of Australia (241944, 339462, 389927, 389875, 389891, 389892, 389938, 443036, 442915, 442981, 496610, 496739, 552485, 552498, 1050208, 1075175). L.Y. is funded by Australian Research Council (Grant number DE200100425). The Minnesota Center for Twin and Family Research (MCTFR) was supported in part by USPHS Grants from the National Institute on Alcohol Abuse and Alcoholism (AA09367 and AA11886) and the National Institute on Drug Abuse (DA05147, DA13240, and DA024417). The Women's Genome Health Study (WGHS) was funded by the National Heart, Lung, and Blood Institute (HL043851 and HL080467) and the National Cancer Institute (CA047988 and UM1CA182913), with support for genotyping provided by Amgen. Data collection in the Finnish Twin Registry has been supported by the Wellcome Trust Sanger Institute, the Broad Institute, ENGAGE-European Network for Genetic and Genomic Epidemiology, FP7-HEALTH-F4-2007, grant agreement number 201413, National Institute of Alcohol Abuse and Alcoholism (grants AA-12502, AA-00145, AA-09203, AA15416, and K02AA018755) and the Academy of Finland (grants 100499, 205585, 118555, 141054, 264146, 308248, 312073 and 336823 to J. Kaprio). TwinsUK is funded by the Wellcome Trust, Medical Research Council, Versus Arthritis, European Union Horizon 2020, Chronic Disease Research Foundation (CDRF), Zoe Ltd and the National Institute for Health Research (NIHR) Clinical Research Network (CRN) and Biomedical Research Centre based at Guy's and St Thomas' NHS Foundation Trust in partnership with King's College London. For NESDA, funding was obtained from the Netherlands Organization for Scientific Research (Geestkracht program grant 10000-1002), the Center for Medical Systems Biology (CSMB, NVVO Genomics), Biobanking and Biomolecular Resources Research Infrastructure (BBMRI-NL), VU University's Institutes for Health and Care Research (EMGO+) and Neuroscience Campus Amsterdam, University Medical Center Groningen, Leiden University Medical Center, National Institutes of Health (NIH, ROI D0042157-01A, MH081802, Grand Opportunity grants 1 RC2 Ml-1089951 and IRC2 MH089995). Part of the genotyping and analyses were funded by the Genetic Association Information Network (GAIN) of the Foundation for the National Institutes of Health. Computing was supported by BiG Grid, the Dutch e-Science Grid, which is financially supported by NWO. Work in the Del Bene lab was supported by the Programme Investissements d'Avenir IHU FOReSIGHT (ANR-18-IAHU-01). C.R. was supported by an EU Horizon 2020 Marie Sklodowska-Curie Action fellowship (H2020-MSCA-IF-2014 #661527). H.S. and K.S. are employees of deCODE Genetics/Amgen. The other authors declare no competing financial interests. TRIAL REGISTRATION NUMBER: N/A.

GWAS of Dizygotic Twinning in an Enlarged Australian Sample of Mothers of DZ Twins.

Female fertility is a complex trait with age-specific changes in spontaneous dizygotic (DZ) twinning and fertility. To elucidate factors regulating female fertility and infertility, we conducted a genome-wide association study (GWAS) on mothers of spontaneous DZ twins (MoDZT) versus controls (3273 cases, 24,009 controls). This is a follow-up study to the Australia/New Zealand (ANZ) component of that previously reported (Mbarek et al., 2016), with a sample size almost twice that of the entire discovery sample meta-analysed in the previous article (and five times the ANZ contribution to that), resulting from newly available additional genotyping and representing a significant increase in power. We compare analyses with and without male controls and show unequivocally that it is better to include male controls who have been screened for recent family history, than to use only female controls. Results from the SNP based GWAS identified four genomewide significant signals, including one novel region, ZFPM1 (Zinc Finger Protein, FOG Family Member 1), on chromosome 16. Previous signals near FSHB (Follicle Stimulating Hormone beta subunit) and SMAD3 (SMAD Family Member 3) were also replicated (Mbarek et al., 2016). We also ran the GWAS with a dominance model that identified a further locus ADRB2 on chr 5. These results have been contributed to the International Twinning Genetics Consortium for inclusion in the next GWAS meta-analysis (Mbarek et al., in press).

Genetic predisposition to porto-sinusoidal vascular disorder: A functional genomic-based, multigenerational family study.

BACKGROUND AND AIMS: Porto-sinusoidal vascular disorder (PSVD) is a group of liver vascular diseases featuring lesions encompassing the portal venules and sinusoids unaccompanied by cirrhosis, irrespective of the presence/absence of portal hypertension. It can occur secondary to coagulation disorders or insult by toxic agents. However, the cause of PSVD remains unknown in most cases. Hereditary cases of PSVD are exceptionally rare, but they are of particular interest and may unveil genetic alterations and molecular mechanisms associated with the disease. APPROACH AND RESULTS: We performed genome sequencing of four patients and two healthy individuals of a large multigenerational Lebanese family with PSVD and identified a heterozygous deleterious variant (c.547C>T, p.R183W) of FCH and double SH3 domains 1 ( FCHSD1 ), an uncharacterized gene, in patients. This variant segregated with the disease, and its pattern of inheritance was suggestive of autosomal dominant with variable expressivity. RNA structural modelling of human FCHSD1 suggests that the C-to-T substitution at position 547, corresponding to FCHSD1R183W , may increase both messenger RNA (mRNA) and protein stability and its interaction with MTOR-associated protein, LST8 homolog, a key protein of the mechanistic target of rapamycin (mTOR pathway). These predictions were substantiated by biochemical analyses, which showed that FCHSD1R183W induced high FCHSD1 mRNA stability, overexpression of FCHSD1 protein, and an increase in mTORC1 activation. This human FCHSD1 variant was introduced into mice through CRISPR/Cas9 genome editing. Nine out of the 15 mice carrying the human FCHSD1R183W variant mimicked the phenotype of human PSVD, including splenomegaly and enlarged portal vein. CONCLUSIONS: Aberrant FCHSD1 structure and function leads to mTOR pathway overactivation and may cause PSVD.

Overcoming Eurocentric bias makes for better science.

To understand disease, scientists are producing comprehensive omics datasets. However, the majority of these are Eurocentric. Recently, the inclusion of patients from Asia and the Middle East in genomic analyses uncovered unique loci linked to COVID-19 severity. This demonstrates that focusing on diversity and underrepresented populations can benefit all.

Analysis of incidental findings in Qatar genome participants reveals novel functional variants in LMNA and DSP.

In order to report clinically actionable incidental findings in genetic testing, the American College of Medical Genetics and Genomics (ACMG) recommended the evaluation of variants in 59 genes associated with highly penetrant mutations. However, there is a lack of epidemiological data on medically actionable rare variants in these genes in Arab populations. We used whole genome sequencing data from 6045 participants from the Qatar Genome Programme and integrated it with phenotypic data collected by the Qatar Biobank. We identified novel putative pathogenic variants in the 59 ACMG genes by filtering previously unrecorded variants based on computational prediction of pathogenicity, variant rarity and segregation evidence. We assessed the phenotypic associations of candidate variants in genes linked to cardiovascular diseases. Finally, we used a zebrafish knockdown and synthetic human mRNA co-injection assay to functionally characterize two of these novel variants. We assessed the zebrafish cardiac function in terms of heart rate, rhythm and hemodynamics, as well as the heart structure. We identified 52 492 novel variants, which have not been reported in global and disease-specific databases. A total of 74 novel variants were selected with potentially pathogenic effect. We prioritized two novel cardiovascular variants, DSP c.1841A > G (p.Asp614Gly) and LMNA c.326 T > G (p.Val109Gly) for functional characterization. Our results showed that both variants resulted in abnormal zebrafish heart rate, rhythm and structure. This study highlights medically actionable variants that are specific to the Middle Eastern Qatari population.

Qatar genome: Insights on genomics from the Middle East.

Despite recent biomedical breakthroughs and large genomic studies growing momentum, the Middle Eastern population, home to over 400 million people, is underrepresented in the human genome variation databases. Here we describe insights from Phase 1 of the Qatar Genome Program with whole genome sequenced 6047 individuals from Qatar. We identified more than 88 million variants of which 24 million are novel and 23 million are singletons. Consistent with the high consanguinity and founder effects in the region, we found that several rare deleterious variants were more common in the Qatari population while others seem to provide protection against diseases and have shaped the genetic architecture of adaptive phenotypes. These results highlight the value of our data as a resource to advance genetic studies in the Arab and neighboring Middle Eastern populations and will significantly boost the current efforts to improve our understanding of global patterns of human variations, human history, and genetic contributions to health and diseases in diverse populations.

The QChip1 knowledgebase and microarray for precision medicine in Qatar.

Risk genes for Mendelian (single-gene) disorders (SGDs) are consistent across populations, but pathogenic risk variants that cause SGDs are typically population-private. The goal was to develop "QChip1," an inexpensive genotyping microarray to comprehensively screen newborns, couples, and patients for SGD risk variants in Qatar, a small nation on the Arabian Peninsula with a high degree of consanguinity. Over 108 variants in 8445 Qatari were identified for inclusion in a genotyping array containing 165,695 probes for 83,542 known and potentially pathogenic variants in 3438 SGDs. QChip1 had a concordance with whole-genome sequencing of 99.1%. Testing of QChip1 with 2707 Qatari genomes identified 32,674 risk variants, an average of 134 pathogenic alleles per Qatari genome. The most common pathogenic variants were those causing homocystinuria (1.12% risk allele frequency), and Stargardt disease (2.07%). The majority (85%) of Qatari SGD pathogenic variants were not present in Western populations such as European American, South Asian American, and African American in New York City and European and Afro-Caribbean in Puerto Rico; and only 50% were observed in a broad collection of data across the Greater Middle East including Kuwait, Iran, and United Arab Emirates. This study demonstrates the feasibility of developing accurate screening tools to identify SGD risk variants in understudied populations, and the need for ancestry-specific SGD screening tools.

One Year of SARS-CoV-2: Genomic Characterization of COVID-19 Outbreak in Qatar.

Qatar, a country with a strong health system and a diverse population consisting mainly of expatriate residents, has experienced two large waves of COVID-19 outbreak. In this study, we report on 2634 SARS-CoV-2 whole-genome sequences from infected patients in Qatar between March-2020 and March-2021, representing 1.5% of all positive cases in this period. Despite the restrictions on international travel, the viruses sampled from the populace of Qatar mirrored nearly the entire global population's genomic diversity with nine predominant viral lineages that were sustained by local transmission chains and the emergence of mutations that are likely to have originated in Qatar. We reported an increased number of mutations and deletions in B.1.1.7 and B.1.351 lineages in a short period. These findings raise the imperative need to continue the ongoing genomic surveillance that has been an integral part of the national response to monitor the SARS-CoV-2 profile and re-emergence in Qatar.

Poking COVID-19: Insights on Genomic Constraints among Immune-Related Genes between Qatari and Italian Populations.

Host genomic information, specifically genomic variations, may characterize susceptibility to disease and identify people with a higher risk of harm, leading to better targeting of care and vaccination. Italy was the epicentre for the spread of COVID-19 in Europe, the first country to go into a national lockdown and has one of the highest COVID-19 associated mortality rates. Qatar, on the other hand has a very low mortality rate. In this study, we compared whole-genome sequencing data of 14398 adults and Qatari-national to 925 Italian individuals. We also included in the comparison whole-exome sequence data from 189 Italian laboratory-confirmed COVID-19 cases. We focused our study on a curated list of 3619 candidate genes involved in innate immunity and host-pathogen interaction. Two population-gene metric scores, the Delta Singleton-Cohort variant score (DSC) and Sum Singleton-Cohort variant score (SSC), were applied to estimate the presence of selective constraints in the Qatari population and in the Italian cohorts. Results based on DSC and SSC metrics demonstrated a different selective pressure on three genes (MUC5AC, ABCA7, FLNA) between Qatari and Italian populations. This study highlighted the genetic differences between Qatari and Italian populations and identified a subset of genes involved in innate immunity and host-pathogen interaction.

Identical twins carry a persistent epigenetic signature of early genome programming.

Monozygotic (MZ) twins and higher-order multiples arise when a zygote splits during pre-implantation stages of development. The mechanisms underpinning this event have remained a mystery. Because MZ twinning rarely runs in families, the leading hypothesis is that it occurs at random. Here, we show that MZ twinning is strongly associated with a stable DNA methylation signature in adult somatic tissues. This signature spans regions near telomeres and centromeres, Polycomb-repressed regions and heterochromatin, genes involved in cell-adhesion, WNT signaling, cell fate, and putative human metastable epialleles. Our study also demonstrates a never-anticipated corollary: because identical twins keep a lifelong molecular signature, we can retrospectively diagnose if a person was conceived as monozygotic twin.

Actionable genomic variants in 6045 participants from the Qatar Genome Program.

In a clinical setting, DNA sequencing can uncover findings unrelated to the purpose of genetic evaluation. The American College of Medical Genetics and Genomics (ACMG) recommends the evaluation and reporting of 59 genes from clinic genomic sequencing. While the prevalence of secondary findings is available from large population studies, these data lack Arab and other Middle Eastern populations. The Qatar Genome Program (QGP) generates whole-genome sequencing (WGS) data and combines it with phenotypic information to create a comprehensive database for studying the Qatari and wider Arab and Middle Eastern populations at the molecular level. This study identified and analyzed medically actionable variants in the 59 ACMG genes using WGS data from 6045 QGP participants. Our results identified a total of 60 pathogenic and likely pathogenic variants in 25 ACMG genes in 141 unique individuals. Overall, 2.3% of the QGP sequenced participants carried a pathogenic or likely pathogenic variant in one of the 59 ACMG genes. We evaluated the QGP phenotype-genotype association of additional nonpathogenic ACMG variants. These variants were found in patients from the Hamad Medical Corporation or reported incidental findings data in Qatar. We found a significant phenotype association for two variants, c.313+3A>C in LDLR, and c.58C>T (p.Gln20*) in the TPM1.

Genetic meta-analysis of twin birth weight shows high genetic correlation with singleton birth weight.

Birth weight (BW) is an important predictor of newborn survival and health and has associations with many adult health outcomes, including cardiometabolic disorders, autoimmune diseases and mental health. On average, twins have a lower BW than singletons as a result of a different pattern of fetal growth and shorter gestational duration. Therefore, investigations into the genetics of BW often exclude data from twins, leading to a reduction in sample size and remaining ambiguities concerning the genetic contribution to BW in twins. In this study, we carried out a genome-wide association meta-analysis of BW in 42 212 twin individuals and found a positive correlation of beta values (Pearson's r = 0.66, 95% confidence interval [CI]: 0.47-0.77) with 150 previously reported genome-wide significant variants for singleton BW. We identified strong positive genetic correlations between BW in twins and numerous anthropometric traits, most notably with BW in singletons (genetic correlation [rg] = 0.92, 95% CI: 0.66-1.18). Genetic correlations of BW in twins with a series of health-related traits closely resembled those previously observed for BW in singletons. Polygenic scores constructed from a genome-wide association study on BW in the UK Biobank demonstrated strong predictive power in a target sample of Dutch twins and singletons. Together, our results indicate that a similar genetic architecture underlies BW in twins and singletons and that future genome-wide studies might benefit from including data from large twin registers.

Rare loss-of-function variants in type I IFN immunity genes are not associated with severe COVID-19.

A recent report found that rare predicted loss-of-function (pLOF) variants across 13 candidate genes in TLR3- and IRF7-dependent type I IFN pathways explain up to 3.5% of severe COVID-19 cases. We performed whole-exome or whole-genome sequencing of 1,864 COVID-19 cases (713 with severe and 1,151 with mild disease) and 15,033 ancestry-matched population controls across 4 independent COVID-19 biobanks. We tested whether rare pLOF variants in these 13 genes were associated with severe COVID-19. We identified only 1 rare pLOF mutation across these genes among 713 cases with severe COVID-19 and observed no enrichment of pLOFs in severe cases compared to population controls or mild COVID-19 cases. We found no evidence of association of rare LOF variants in the 13 candidate genes with severe COVID-19 outcomes.

Loss of C2orf69 defines a fatal autoinflammatory syndrome in humans and zebrafish that evokes a glycogen-storage-associated mitochondriopathy.

Human C2orf69 is an evolutionarily conserved gene whose function is unknown. Here, we report eight unrelated families from which 20 children presented with a fatal syndrome consisting of severe autoinflammation and progredient leukoencephalopathy with recurrent seizures; 12 of these subjects, whose DNA was available, segregated homozygous loss-of-function C2orf69 variants. C2ORF69 bears homology to esterase enzymes, and orthologs can be found in most eukaryotic genomes, including that of unicellular phytoplankton. We found that endogenous C2ORF69 (1) is loosely bound to mitochondria, (2) affects mitochondrial membrane potential and oxidative respiration in cultured neurons, and (3) controls the levels of the glycogen branching enzyme 1 (GBE1) consistent with a glycogen-storage-associated mitochondriopathy. We show that CRISPR-Cas9-mediated inactivation of zebrafish C2orf69 results in lethality by 8 months of age due to spontaneous epileptic seizures, which is preceded by persistent brain inflammation. Collectively, our results delineate an autoinflammatory Mendelian disorder of C2orf69 deficiency that disrupts the development/homeostasis of the immune and central nervous systems.

Genome-wide association study of panic disorder reveals genetic overlap with neuroticism and depression.

Panic disorder (PD) has a lifetime prevalence of 2-4% and heritability estimates of 40%. The contributory genetic variants remain largely unknown, with few and inconsistent loci having been reported. The present report describes the largest genome-wide association study (GWAS) of PD to date comprising genome-wide genotype data of 2248 clinically well-characterized PD patients and 7992 ethnically matched controls. The samples originated from four European countries (Denmark, Estonia, Germany, and Sweden). Standard GWAS quality control procedures were conducted on each individual dataset, and imputation was performed using the 1000 Genomes Project reference panel. A meta-analysis was then performed using the Ricopili pipeline. No genome-wide significant locus was identified. Leave-one-out analyses generated highly significant polygenic risk scores (PRS) (explained variance of up to 2.6%). Linkage disequilibrium (LD) score regression analysis of the GWAS data showed that the estimated heritability for PD was 28.0-34.2%. After correction for multiple testing, a significant genetic correlation was found between PD and major depressive disorder, depressive symptoms, and neuroticism. A total of 255 single-nucleotide polymorphisms (SNPs) with p < 1 × 10-4 were followed up in an independent sample of 2408 PD patients and 228,470 controls from Denmark, Iceland and the Netherlands. In the combined analysis, SNP rs144783209 showed the strongest association with PD (pcomb = 3.10 × 10-7). Sign tests revealed a significant enrichment of SNPs with a discovery p-value of <0.0001 in the combined follow up cohort (p = 0.048). The present integrative analysis represents a major step towards the elucidation of the genetic susceptibility to PD.